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Abstract. This is the third paper of a series of papers where we explore the evolution of iron-rich ejecta from quark-novae. In 
the first paper, we explored the case where a quark-nova ejecta forms a degenerate shell, supported by the star's magnetic field, 
with applications to SGRs. In the second paper we considered the case where the ejecta would have sufficient angular momen- 
tum to form a degenerate Keplerian torus and applied such a system to two AXPs, namely 1E2259+586 and 4U0142+615. Here 
we explore the late evolution of the degenerate torus and find that it can remain unchanged for ~ 10 6 years before it becomes 
non-degenerate. This transition from a degenerate torus (accretion dominated) to a non-degenerate disk (no accretion), occurs 
about 10 6 years following the quark-nova, and exhibits features that are reminiscent of observed properties of RRATs. Using 
this model, we can account for the duration of both the radio bursts and the quiet phase, as well as the observed radio flux 
from RRATs. The unique on and off activity of the radio pulsar PSR B 1 93 1 +24 is similar to that of "old RRATs" in our model. 
For old RRATs, in our model, the spin-down rate during the radio loud phase is about a factor 1.6 larger than when it is quiet, 
remarkably similar to what has been measured for PSR B 193 1+24. We discuss a connection between XDINs and RRATs and 
argue that some XDINs may be "dead RRATs" that have already consumed their non-degenerate disk. 

Key words, accretion, accretion disks - (stars:) pulsars: general - dense matter - X-rays: bursts - Elementary particles 



1. Introduction 



|McLaughlin et al. (2006) have recently reported the detection 
of eleven "Rotation RAdio Transients", or "RRATs", charac- 
terised by repeated, irregular radio bursts, with burst dura- 
tions of 2-30 ms, and intervals between bursts of ~ 4 min 
to ~ 3 hr. The RRATs are concentrated at low Galactic lat- 
itudes, with distances implied by their dispersion measures 
of ~ 2-7 kpc. For ten of the eleven RRATs discovered by 
|McLaughlin et al. (2006)| an analysis of the spacings between 
repeat bursts reveals an underlying spin period, P, and also in 
three cases, a spin period derivative, P. The observed periods 
fall in the range 0.4 s < P < 7 s, which generally overlap 
with those seen for the radio pulsar population. Since August 
2003, all the sources have been reobserved at least nine times 
at intervals of between one and six months (i.e. they show spo- 
radic radio bursting for years). All have shown multiple bursts, 
with between 4 and 229 events detected in total from each ob- 
ject (see Table 1 in McLaughlin et al. (2006)1. For the three 



The discovery of the first X-ray counterpart to a RRAT 
was reported by |Reynolds et al. (2006)| for RRAT J1819-1458. 
The X-ray spectrum detected by Chandra and XMM-Newton is 
thermal, fit well by a blackbody with kT ~ 0.1 keV with a pos- 
sible excess at high energies. Combined with the high inferred 
surface magnetic field strength, long spin period, and lack of 
persistent radio emission this led to a comparison with the pop- 



ulation of magnetars (e.g. R eynolds et al. (2006)| l. These com- 
parative studies showed the obvious differences between mag- 
netars and RRATs: (i) the magnetar birth rate is well below that 



RRATs with values measured for both P and P, a characteristic 
age, t c , and a dipole surface magnetic field, B, can be inferred 
and are listed in Table [1] The total number of these objects 
is a few times that previously estimated for all radio pulsars 



estimated for the RRAT population ( Popov etal. (2006)); (ii) 
X-ray properties of RRATs seem distinct from those seen for 
magnetars: RRATs seem much colder and less luminous than 
the magnetars, and apparently lacks the hard X-ray tail seen 
for these sources (e.g. |Gaensler et al. (2007)l|Rea (2007)1 ; (iii) 
X-ray temperature of RRATs is around 0.1-0.2 keV below the 
average 0.5 keV associated with magnetars. 

The lack of persistent radio emission and the long spin 
period of at least one RRAT has raised the possibility of a 
link between XDINs (Haberl 2004^ and RRATs. Furthermore, 



Popov et al. (2006) show that the inferred birthrate of RRATs 



(McLaughlin et al. (2006) i. 



is consistent with that of XDINs but not with magnetars. The 
XDINs are slightly cooler (kT ~ 0.04 - 0. 1 keV) and less lumi- 
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1 No radio emission of any kind has been reported from XDINs; 
recent observations (Bradley 2006) show no RRAT-like radio bursts 
toward RX J0720.4— 3125 (or toward magnetars). 
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nous (L x ~ 10 31 - 10 32 ergs s" 1 ) than J1819-1458. However, 
the measured period derivatives of two XDINs (RBS 1223 and 
RX J0720.4-3125; Kaplan & van Kerkwijk 2005a, 2005b), 
and the detection of possible proton cyclotron lines in their 
spectra (van Kerkwijk 2004), imply magnetic field strengths 
similar to those of J1819-1458. The immediate questions are 
(i) why only one RRAT shows X-rays as the XDINs do, 
and (ii) why XDINs do not show any radio emission (latest 
searches for pulsed and bursty radio emission from XDINs led 
to no detection despite of their proximity compared to RRATs; 
|Kondratiev et al. (2007)| ). 

RRAT J18 19— 1458 properties in the X-ray show similari- 
ties to those of radio pulsars with ages around 100 kyr. For ex- 
ample, PSR J0538+2817 is 30 kyr old and has kT m = 160 eV, 
while PSR B0656+14 is 110 kyr old and has kT^ = 70 eV 



(see Reynolds et al. (2006) for details). However, the inferred 
surface magnetic field strength of RRATs is at least an order 
of magnitude greater than the radio sources. Two radio pul- 
sars with comparable magnetic fields that have been detected in 
X-rays are PSRs J17 18-37 18 (|Kaspi and McLaughlin (2005)] l 
and Jl 119-6127 (Gon zalez et al. (2005) I. These sources show 
temperatures (kT ~ 150 - 200 eV) and luminosities (~ 
10 32 - 10 33 ergs s 1 ) comparable to that of RRAT J1819-1458, 
although both sources are probably much younger (35 and 
1.7 kyr, respectively) and, in contrast to RRAT J1819-1456, 
have X-ray luminosities less than their spin-down luminosities. 



1.1. Literature explanations 

The models proposed so far in the literature can be classi- 
fied into the following categories: (i) Extreme pulses from 
distant pulsars ( |Weltevrede et al.(2006)| l, similar to the pulses 
seen from the nearby pulsar B0656+14. (ii) Re-activated ra- 
dio pulsars ( Zhang et al. (2007)} where the RRATs are pul- 
sars that are no longer active, but for which a temporary "star 
spot" with multipole field components emerges above the sur- 
face. This magnetic field component could temporarily reac- 
tivate the radio beaming mechanism, producing the observed 
bursts, (iii) Nulling pulsars viewed from the opposite direc- 
tion (Zh ang et al. (2007)] l. In this model, RRATs are normal 
pulsars with their magnetic poles not aligned favorably for de- 
tection, but undergo an magnetic reversal and so occasionally 
produce emission that can be observed, (iv) Sporadic accre- 
tion ( |Cordes and Shannon (2006)] |Li (2006)| l. Here the RRAT 
mechanism might be produced by interaction of the neutron 
star with an equatorial fallback disk or with orbiting circum- 
polar debris. Accretion from a disk should quench the radio 
emission mechanism, but sporadic drops in the accretion rate 
could allow the radio beam to turn on for a fraction of a sec- 
ond, producing the RRAT phenomenon (Li (2006)). 

In this paper, the third of a series, we present an alternative 
scenario that incorporates SGRs, AXPs, XDINs, and RRATs 
into one family of compact objects. In our model, SGRs, AXPs, 
XDINs, and RRATs are all strange quark stars that differ, for 
the most part, only by age. 

This paper is organized as follows: Sect. 2 describes our 
strange quark star model. Sect. 3 discusses the implications of 



the evolution of debris material left over from the birth of the 
strange quark star. In Sect. 4 we show how after roughly one 
million years the debris can become responsible for RRAT be- 
haviour. Then, in Sect. 5 we describe observations and make 
predictions using our model. We discuss some further implica- 
tions of our model in Sect. 5, and lastly we conclude in Sect. 6. 

2. Our model 



In the quark-nova (QN) picture (|Ouyed et al. (2002) 



Kera nen et al. (2005){ hereafter ODD and KOJ respectively) 
the core of a neutron star, that undergoes the phase transition to 
the quark phase, shrinks in a spherically symmetric fashion to 
a stable, more compact strange matter configuration faster than 
the overlaying material (the neutron-rich hadronic envelope) 
can respond, leading to an effective core collapse. The core 
of the neutron star is a few kilometers in radius initially, 
but shrinks to 1-2 km in a collapse time of about 0.1 ms 



(Lugones et al. (1994) i. The gravitational potential energy 
released (plus latent heat of phase transition) during this event 
is converted partly into internal energy and partly into an 
outward propagating shock wave that imparts kinetic energy to 
the overlying material. 



As described in a series of papers ( Ouyed etal. (2007a) 



|Ouyed et al. (2007b){ hereafter referred to as OLNI and 
OLNII), during a quark-nova the degenerate crust of a neu- 
tron star is blown off, leaving behind a strange quark star (QS) 
surrounded by left over, highly-metallic degenerate matter. In 
OLNI we discussed the case where the ejected crust had insuffi- 
cient angular momentum to escape the QS's gravitational pull, 
and so would either balance with the QS's magnetic field and 
form a co-rotating shell, or fall back entirely. Then in OLNII 
the case where the QS was born with a sufficient spin-period 
to impart the ejected crust with enough angular momentum to 
form a degenerate torus. In this paper, we explore the result 
of this degenerate torus after enough time has passed for it to 
expand to densities where it is non-degenerate. 

2.1. Emission: Vortex Expulsion vs. Accretion 

In our model the QS is in the ground Color-Flavor Locked 
(CFL) phase and so behaves as a type II superconductor, 
wherein a rotationally-induced lattice forms inside the star. As 
the star spins down the magnetic field, which is confined to ex- 
ist only within the vortices, is also expelled and the subsequent 
magnetic reconnection leads to the production of X-rays. The 
luminosity from vortex expulsion is given in OLNI to be, 



Lx,v - 2.01 x 10 34 erg S-V0.1PI11 • 



(1) 



Here the subscript v stands for "vortex" and 77X.0.1 is the ef- 
ficiency parameter, in units of 0.1, inherent in the conversion 
from magnetic energy to observed radiation. The spin-down 
rate, P is given in units of 10~ n s s -1 . 

In this paper the degenerate material ejected during the 
quark-nova is imparted with sufficient angular momentum to 
form a degenerate torus, in which case accretion from this torus 
can result in emission that outshines the emission due to vortex 
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expulsion. The condition on the initial spin-period of the QS at 
birth is then an upper limit (from OLNII), 



Po < Po., 



2.5 ms 



R 3/2 ff 9/2 
D 0,15 n QS,10 

m'J\MV A 



(2) 



0-7'"QS,1.4 

where Aqs.io is the radius of the QS in units of 10 km, Z?o,i5 
is the initial magnetic field strength in units of 10 15 G, and the 
initial mass of the torus is ;«o._7 in units of 10~ 7 solar masses. 
From here on, we parametrize the actual period as being some 
fraction, ao.3, of the maximum period in units of 0.3; Pq = 

aPo,max- 

The accretion luminosity from the type of torus considered 
in this paper is given by equation (13) in OLNII, 



ix,acc. =< 1.7 xlO 33 ergs" 1 



^3.3 QS 



(3) 



Here 770.1 is the efficiency of conversion of accreted material 
into X-ray emission in units of 0.1, fl t ,i5 is the radial distance of 
the torus in units of 15 km, yi/3.3 is the mean molecular weight 
per electron in units of 3.3 (the quiescent phase value), and 
Mqs,ia is the mass of the QS in units of 1 .4 solar masses. 

2.2. Degenerate Torus properties 

The radius of the torus can be found from the initial spin-period 
of the QS as was done in OLNII, 

^t=15kma- 8 /X 3 s M^ 14 , (4) 

where the initial period is in units of milliseconds. 

By assuming a constant accretion rate, the evolution of the 
torus mass can be determined, and the time needed for the torus 
to reach a density where it becomes non-degenerate is, 

W ° iA ,nS m °'- 7M QS,1.4^3.3 

T t 3.4 x 10 yrs 7— —, , (5) 

with m t being the torus accretion rate given by eq(10) in OLNII. 
The mass and density of the torus over time are plotted in Fig. 

ru 

While the torus is above degenerate densities, diffusion by 
the QS's magnetic field into the inner walls of the torus (lead- 
ing to bursting accretion events; see OLNII) increases the inner 
radius. We can estimate the change inner radius of the torus, 
R'm = Rt + A/? t (see Appendix A), by using the typical time 
needed for magnetic field diffusion into the torus, tb, and the 
typical depth into the torus at which shear stresses become 
great enough such that accretion can proceed, Ar w . Both of 
these parameters were estimated in OLNII (eqns. 17 & 5). We 
find that the change in inner radius over the span of the lifetime 
of the torus is negligible. 

While the magnetic field slowly consumes the torus' inner 
edge, its outer edge moves outward (due to viscosity from par- 
ticle collisions within a degenerate ideal gas) at a rate given by 
equation (A. 7) in OLNII. After the time where the torus ex- 
pands to densities at which it becomes non-degenerate, r t , the 
torus will have extended to an outer radius of, 



flout - 207 km 



in 12 M 2 u 3 



3/2^7/4 



(6) 
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Fig. 1. The mass and density of the torus, assuming a temper- 
ature of 1 keV. T t and p n( j are the time and density at which the 
torus becomes non-degenerate. 



3. Transition to RRAT phase 

The evolution of the density of the torus (during its degenerate 
phase) is calculated using eq(A.8) in OLNII, and is shown in 
the lower panel in Fig. Q] At time r t , the torus density drops 
rapidly and eventually reaches densities that are below de- 
generacy. For temperatures estimated in the torus (eq. 13 in 
ONLII; r eV - 857/n.ii?t,i5/A*QS,i.4), tne density at which the 
torus makes the transition from degeneracy to non-degeneracy 
is, p nd - 1 1.4 gm cm- 3 VqM^KIia- 

Upon the transition from degeneracy to non-degeneracy, 
the new viscosity becomes v = 6.0 x lO -4 !^^ 2 cm 2 s -1 
(OLNII eq. A. 5). The corresponding accretion rate then is 
m ~ 0.0057 g s _1 PndHdT^ 2 which is completely negligible 



(Fra nk et al. (2002)] ); Hi is the disk thickness in units of cen- 
timeters. Thus, immediately following this transition the X-ray 
luminosity of the system is no longer dominated by accretion 
but rather by vortex (magnetic flux) expulsion from the star. 
The consequence of this would be a decrease in X-ray luminos- 
ity by a few orders of magnitude in a relatively short amount of 
time (see left panel in Figure 2). 

3. 1 . Stellar properties 

The evolution of the star's period and period deriva- 
tive are defined by the physics of vortex expulsion (see 



Niebergal et al. (2006)) 



Pqs (t) = Pqs,o 



1 + — 

Tv 



(7) 
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and 



Pqs U) 



where 



r QS,0 
3t v 



1 + 



T V = 840 S BQs,o,15^QS,0,m S ^QS,1.4^Q S , 



10' 



(8) 



(9) 



is the characteristic age due to vortex expulsion. The other two 
important relations in our model are (see OLNI and OLNII) 



B = ^3kPP 
PB 2 = P B 2 Q . 



(10) 



where k = 8.8 x 10 G s . The former equation describes 
magnetic field decay while the latter allows us to link the initial 
and current conditions. 

By making use of the above equations and equation @, the 
star's period, period derivative, and magnetic field can be deter- 
mined at the time, r t , when the transition to a non-degenerate 
torus is made (with r t >> r v ), 



2.3 s 



^q,3.3 



R 4/3 ,,1/2 „10/3 
D 0,15' ra QS,1.4 n QS,10 



7.3 x 10~ 14 s s- 1 P 



-l 



% 1^1,1 5 



Combining 

#0,15 - 



equation 

2 1 f p35/38 pl3/38 



■r B 2/3 M 2 R 32 ' 3 

^q,3.3 D 0,15 /K QS,1.4"QS,10 

above with eq.dTOb 
with f trans . 



trans. ,3 trans. ,-13 



we 
and 



(11) 



find 



units of 3 s and 10 



-13 



r l/3 -6/13 8/3S-39/76 

■ •■' / (.)s.i.-i / Au "o. s > 

Equations (O and ( fTTT i give us R t ^ and Pq 4 



respectively, while 
is of order unity. 

s- 

With P Q , m s, Bo.15 

and eq© together, we solve for mo, -7. 

Table 2 shows birth parameters of the 3 RRATs studied here 
for a = 0.3. We find an initial magnetic field 10 14 < Bo G < 
10 16 and sub-millisecond birth periods. Also listed in Table 2 
are the torus' initial inner radius, R t , and initial mass, mo which 
implies 1O" 8 M < m < 10~ 5 M Q and 10 < R t (km) < 35. 

3.2. Non-degenerate disk properties 

The non-degenerate disk's thickness at a radial distance r 
is given by Ha = v 2 h /g, where g = GMqsH^/ 'r 3 is the 



effective gravity at r , and v t h = yjkT/(ii q mn) 
10 7 cm s r^T^y, is the disk thermal speed. Thus, atR 
find 



7-1/2 ff 3/2 
„ „ , „4 d.keV in.100 

ff d ~ 3.9 x 10 4 cm —j- -fr- 
izz tji/j. 

^q,3.3-™QS,1.4 



1.7 X 

n, we 
(12) 



while the disk's thickness is a few hundred meters at its outer 
edge, R out . However, as we show below, the disk quickly cools 
reducing its temperature and thus its thickness until it solidifies. 
The density would continue dropping below p n( j until it reaches 
the density of normal iron ppe ~ 7.86 g cm -3 . 
The thermal evolution of the disk is given by, 

W - ^- - ~. — L X,v - i^BB.t > 

ot An 



where = (InHd/Rd) is the solid angle extended by the disk 
with thickness In the above, C v is the torus heat capacity 
while LfiB.t — 27r/?g Ut crr 4 is the torus blackbody cooling. The 
disk cools, while Ha decreases, yielding smaller Qd thus further 
cooling. However, the decrease in Ha stops when the disk gas 
condenses which occurs at a temperature of rpe,cond. — 0.265 
eV (CRC tables), since the disk composition is dominated by 
iron (see OLNI and OLNII). 

To find the disk's thickness at the condensation point 
we impose conservation of surface density so that HaP r a - 
Hd, con d.PFe ; we assume the surface density constant during the 
gas phase since the gas phase is short-lived so that 



Ha,: 



R 



3/2 



cond. 



0.12 cm 



,100 



1/2 M l/2 
^q,3-3 QS.1.4 



(14) 



The solid disk is extremely thin varying from a millimeter (pa- 
per thin) at up to a few centimeters at its outer edge. 

Using Qd.cond. = (27r//d,cond./^in), we estimate the equilib- 
rium temperature of the condensed disk to be 



T d ~ 18.6 meV- 



R 



pW 1/4 

a -i3''x,o.i 



R 



1/2 



y 1/8 M 1/8 

00(Fq,3.3 J "QS,1.4 



(15) 



where R out is given in units of 1000 km. For the 3 RRATs in the 
order listed listed in Tables we get, ~ 6 meV, ~ 86 meV and 
~ 5.7 meV, respectively. At temperatures below 0.1 eV, con- 
densed iron is in the form of ferrite, or a-iron, a body-centered 
cubic structure. 



4. RRAT-like sporadic behavior 

The inner disk will be slowly penetrated by the QS's magnetic 
field on timescales determined by the induction equation, 



8B 

at 



c 

4-7TCT 



V 2 B . 



(16) 



Here, <x = n e th e z /l e /(m e v lms ) and A e = l/{n e ^<r T ), where n ejth 
is the number density of thermal electrons in the disk, cr T is the 
Thompson scattering cross-section, and the root-mean-square 
electron velocity is v rms = v t h. Therefore, the time needed for 
the magnetic field to penetrate to a depth of the order of Z/d.cond., 



tb ~ 1.3 hrs 



R 3 



in, 100 



T l/2 1/2 ' 
1 d,meV^q,3.3 



(17) 



Once penetrated the inner disk is accreted, on free-fall 
timescales, onto the star along the magnetic field lines with 



D R 3/2 

Rm _ r, in.100 

kcc = 2.3 ms 



Vff 



M 1/2 

JH QS,1.4 



(18) 



(13) where Vff = ^JtHMq^JR^ is the free-fall velocity. 
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Star 






Radio 




Source 


P(s)' 


P(W~ is s s-') 1 


L x (ergs- 1 ) 2 


ton (ms) 1 


'off (hr) 1 


£*adio (erg s" 1 ) 3 


RRAT J 13 17-5759 


2.6 


0.126 


< 7.5 x 10 32 


10 


0.22 


1.4 x 10 31 


RRAT J1819— 1458 


4.3 


5.76 


3.3 x 10 33 


3 


0.057 


5.6 x 10 31 


RRATJ1913+3333 


0.92 


0.0787 


< 9.4 x 10 34 


2 


0.21 


2.5 x 10 31 



'FromjMcLaughlin et al. (2006) 2 The observed X-ray luminosity (Reynolds et al. 2006); 3 The observed radio luminosity assuming a 

bandwidth of 1 GHz JRea (2007)) 



Table 2. Birth parameters predicted by our model (for a = 0.3). 







Star 




Torus 




Source 


Bo(10 13 G) 


Po (ms) 


P (lO -6 s s" 1 ) 


R t (km) 


m (10-'M Q ) 


RRAT 11317-5759 


0.91 


0.27 


1.2 


13.4 


3.3 


RRAT Jl 8 19-1458 


5.4 


0.97 


12 


31.5 


20 


RRAT J1913+3333 


0.30 


0.20 


0.18 


10.8 


0.53 



Table 3. RRATs era parameters predicted by our model (for a = 0.3). 





Star 






Disk 






Radio 




Source 


Lx,v (erg s~') 


Tbb (eV) 


R m (km) 


tf out (km) T A 


(meV) 


ton = 4cc. (ms) 


tog = t b (hr) 


^Radio (erg s -1 ) 


RRAT 11317-5759 


3.19 x 10 28 


14 


96 


1082 


6.0 


2.2 


0.46 


3.3 x 10 31 


RRAT J1819-1458 


6.7 x 10 31 


100 


33 


203 


86 


0.44 


0.005 


2.3 x 10 30 


RRATJ1913+3333 


1.25 x 10 28 


11 


132 


820 


5.7 


3.5 


1.2 


7.1 x 10 31 



Ly from eqdlb. Tnn from eqJ20t. R<„ from ea iA.H , R nu , from eqJ6l>. T,\ from eqdl5b. t ■„■,■ from eqdl8t. tr from eqd!7t, L, af ij n from eqj 19b with 

A cas , = 10 2 . 



Table 4. (Torus and RRATs) Era duration predicted by our model (for a — 0.3). 



Source 


Ttoms (yrs) 


Trial (yrs) 


RRAT J1317-5759 


2.2 x 10 6 


7.9 x 10 5 


RRAT J1819-1458 


7.9 x 10 4 


3.2 x 10 3 


RRAT J1913+3333 


1.3 x 10 s 


4.1 x 10 5 



T,n,,„ from eq.d5b; t,,,„ from eq.d24ll. 



4.1. Radio emission mechanism and fluxes 

Radio emission in our model is triggered by the accretion of the 
inner ring of mass Amj ~ 2^/?; n //^ cond pFe- Each iron nucleon 
accreted leads to the generation of a GeV photon thus triger- 
ing a photon-pair cascade amplifying the accretion energy by a 
factor A casc .. 

The luminosity in radio is then L ra dio = 
A casc .(0.1Am c 2 )/f acc ., or, 

R 5/2 

iradio ~ A cas , 3.7 x 10 29 erg s" 1 i -^ T ~ . (19) 

fl q ,3.3M Qsl4 

As shown in Table 3, the observed radio fluxes can be ac- 
counted for consistently for all 3 candidates if an amplification 
factor of the order of 10 2 is assumed (a multiplicity number for 
secondary pairs can be as high as 10 5 ; e.g. Melrose 1995 and 
references therein). We expect the medium to be transparent 
to radio since it would most likely be emitted much above the 
polar region. We assume all radio bursts originate at the same 
longitude on the quark star. This guarantees that the bursts are 
modulated at the spin period, consistent with the assumption 
required to derive P and P from the observations of RRATs. 



Table 3 lists parameters derived from our model as com- 
pared to observations which shows some encouraging agree- 
ment with measured radio burst timescale and quiet phases. 
The case of RRAT J18 19— 1458 shows smaller bursting time 
and quiet time, as well as smaller radio fluxes probably because 
our model somehow predicts a smaller R m than the actual one. 
For example, using ~ 80 km instead of 33 km (i.e. a ~ 0.2 
instead of a — 0.3) we find much better agreement. 

5. Model Predictions 

5.1. The two blackbodies 

One of the key aspects of the transition from degenerate torus to 
a non-degenerate disk is the shutting-off of accretion thus elim- 
inating one of the two BBs inherent to the accretion era (see 
OLNII). RRATs should show one dominant blackbody emis- 
sion at a temperature given by 

Tbb = (L x ,v/(47rfl 2 s (r)) I/4 ~ 0.04 keV/^ 2 , (20) 

where P is in units of 10~ 13 s s _1 , and a dimmer one at rep- 
resenting the remnant iron disk. In particular, for RRAT J 1 8 1 9- 
1458 for which P_ n = 5.76 we predict r BB ~ 0.1 keV and 
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another blackbody in infrared at ~ 940 K. The blackbody tem- 
perature for the other two candidates are listed in Table 3 which 
shows a BB temperature in the ultra-violet region and disk tem- 
peratures in the Infrared, 50-100 K. 

5.2. Absorption lines 

Since the disk is cooler than the star (7d < PbbX it will act as 
an absorber of the blackbody emission. The column density in 
the inner part of the disk is 



Nd,Fe = 



#d,cond.PFe 

56mn 



10 22 gm cm~ 2 R V2 



100 



(21) 



high enough to lead to absorption. The X-ray spectrum of 
RRAT J18 19— 1458 showed a single blackbody component did 
not fit the data because of the presence of two strong absorp- 
tion features around 0.5 and 1 keV (Rea (2007)). In our model, 
these may be explained as absorption lines due the disk. 

5.3. Increase in spin-down rate when radio is on 

During the quiet period the star's spin-down rate is simply 
given by P = P/(3f age ) where f age is the source age. However, it 
is natural to expect spin-down increase during the on era since 
the polar wind (secondary pairs) act as a torque against the 
star's rotation. 

The angular momentum per unit mass lost at the light cylin- 
der is c 2 /Qqs, which enhances the spin-down rate of the quark 
star during radio bursts to, using (Amd/f acc .)c 2 /QQs - /qs^qs, 
Iqs is the star's moment of inertia, 



pl R 5/2 



AP 6.7 x 10~ w s s" 



in,100 



M QS ,i. 4 /f QSil0 



(22) 



Compared to P t rans. we get 



AP 



9.2 X 10 



p3r>5/2 R 2/3 
-4 r in,100 D 0,15 



Po, ms R 



13/3 
t,15 



(23) 



As can be seen from the equation above the increase in spin- 
down rate when radio is on is negligible is negligible. However, 
as we show below, as the RRAT ages the increase in spin-down 
rate during the on period is noticeable. 

5.4. Old RRATs: the case ofPSR B1931+24 

A simple estimate of the lifetime of a RRAT can be found by 
integrating equation dt/dr =2 tb /Hd,cond., from P; n to R out , to 
get 



Tn-at ~ 1.6 x 10 yrs 



R V2 
out,1000 

d,meV 



(24) 



The late stages of the RRAT (j rral ~ 10 6 yrs) era, correspond to 
the era when the magnetic field is consuming the last outermost 
remnants of the disk; i.e. when the disk radius is close to R out . 

When the radio is off (no accretion) the star spins-down 
via vortex expulsion with a period derivative given by P — 



P/OT mt ). Using equations (l24t and ( fT3T > with P; n ~ R out , we 
get 



off 



2.3 x 10" 14 s s-'P^'^R 



lp4/3p -43/12 



out, 1000 



(25) 



The increase in spin-down rate during the on period, is found 
by replacing R m ~ R out in eq.(l22l. so that 



3p5/2 



AP 



2.1 x 10" 14 s s- 1 



P 3 R 



out, 1000 



M QS ,1.4# QS io 



(26) 



The spin-down rate during the on era is thus P on = P D ff + AP. 

The observed periodicity of the radio-on 
and radio-off recurrence of PSR B 193 1+24 



(Kramer et al. (2006)Kramer, Lyne, O'Brien, Jordan, and Lorimen 
is difficult to explain in any scenario 
considering an isolated pulsar, and as 



Kramer et al. (2006)Kramer, Lyne, O'Brien, Jordan, and Lorimer 



pointed out, the short shut-off time of less than 10 seconds, 
rules out possible scenarios like precession of the neutron star. 
Furthermore, the facts that the off periods of PSR B 193 1+24 
are five orders of magnitude longer than typical nulling 
periods, that the activity pattern is quasi-periodic and that not 
a single null has been observed during on periods strongly 
argues against the nulling scenario. It appears however that 
magnetospheric conditions are sufficient to explain the change 
in the neutron star torque, but it is not clear what determines 
the ~ 30-40 days periodicity or what could be responsible for 
changing the plasma flow in the magnetosphere, in particular 
in this quasi-periodic fashion. 

Assuming PSR B1931+24 (P = 0.813 s) is an old RRAT 
that is currently consuming the outermost part of its disk, our 
model implies a corresponding spin-down frequency during the 
off period of v of f = -v 2 P ff - -2.6 x 10~ 14 Hz s -1 , while 
during the on period the increase in spin-down frequency is 
Av =* -1.7x 10" 14 Hz s . It means an increase of about 65% in 
spin-down rate during the on period so that v on ^ -4.3 x 10" 14 
Hz s . Our numbers are remarkably similar to what has been 
measured for PSR B 193 1+24. For this source the observed 
spin-down rate the pulsar rotation slows down 50% faster when 
it is on than when it is off with v jf = - 1 .08 x 1 0~ 14 Hz s~ 1 while 
v on = -1.63 x 10~ 14 Hz s _1 . 

For PSR B 193 1+24 the quiet period last for ~ 30-40 days 
while the active period is of the order of ~ 5 days. In our 
model, the timescales for the bursting and quiet phases in old 
RRATs are f acc ~ 70 ms and tb ~ 50 days. While the off period 
is again remarkably similar to the observed one, the observed 
on period of ~ 5 days is difficult to reconcile with our model 
which predicts ~ 70 ms. It is worth noticing however that 



Kramer et al. (2006)Kramer, Lyne, O'Brien, Jordan, and Lorimer 
have been able to observe one switch from on to off that oc- 
curred within less than 10 seconds (although switches between 
the on and off states are rare events). This number is rather 
closer to f acc in our model than to the radio-on timescale. It 
might be possible that another mechanism is at play delaying 
the accretion onto the pole. Finally, the corresponding radio 
flux in our model is estimated to be L,- a dio ~ 10 34 



erg s 
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5.5. XDINs and "dead" RRATs in our model 

Some XDINs might be direct descendants of QSs with shells 
(SGRs in our model; see Appendix B). These would have 
evolved along the "vortex" band since their birth. However, as 
illustrated in Figure [2] some XDINs might be dead RRATs. 
That is, they fell into the "vortex" band following torus 
consumption and eventually evolved along the band. The 
main difference between SGR-descendents XDINs and AXP- 
descendent XDINs (dead RRATs) is the possibility of remnant 
disk surrounding the AXP-descendent XDINs. We thus expect 
some of the XDINs to share some common properties with 
RRATs such as optical (or Infrared) excess. Broad absorption 
lines, similar to those seen for RRAT J1819-1458, have been 
observed for six out of seven XDINs (van Kerkwijk & Kaplan 
2007; Haberl 2007). We have already argued in OLNII that 
these absorption lines are caused by absorption from an old, 
cold disk (see also §5.2). 

At first glance the birth rates of RRATs as descendants of 
AXPs in our model, appears to be too low to explain the high 
inferred population (McLaughlin et al. 2006) of a few times 
(~ 5) that of radio pulsars. Taking the RRAT birth rate to be 
that of AXPs, i.e. ~ 1 /(500 yrs) (Leahy&Ouyed 2007 and ref- 
erences therein), implies a birth rate of ~ 1/10 of that of radio 
pulsars. In order to have a RRAT population ~ 5 times the ra- 
dio pulsar population, the RRAT lifetime must be ~ 50 times 
longer than that of a radio pulsar (~ 10 6 yrs). In our model, the 
RRAT lifetime estimate, given by eq.(l24b. is of order a few mil- 
lion years. However, this is strongly dependent on R out . An R out 
larger by a factor of 4 (caused by different effects; see eq.(6)) 
gives a large enough RRAT lifetime, i.e. a low enough birth 
rate. 

6. Conclusion 

Here we have presented a model where RRATs fit in natu- 
rally as descendants of AXPs with degenerate iron-rich disks, 
with the transition occurring when the disk becomes non- 
degenerate. These AXPs in our model are the quark star rem- 
nants of quark-novae, surrounded by the ejected former neu- 
tron star's crust material. The evolution of quark stars as they 
spin down is summarized in the Lx-P diagram in Figure 2 (see 
appendix B). Although this model is speculative, it can ex- 
plain many features of SGRs and AXPs (OLNI and OLNII) 
and RRATs. 
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Appendix A: Evolution of torus inner radius 

During the degenerate phase, the inner walls of the torus are 
carved out by magnetic penetration (leading to bursting accre- 
tion events; see OLNII) increasing the inner radius from R t to 
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P_u ts/s] 

Fig. 2. These figures depict the important regimes of SGRs, AXPs, RRATs, and XDINs in our model (see app. B), the left panel 
being an evolutionary illustration and the right panel contains observational data. There are two important transitions to note 
for AXPs (torus-harbouring quark stars) in our model, the first occurs when the star spins-down to P « 10~ n s/s. At this point 
the dominant luminosity mechanism changes from being magnetic flux (vortex) expulsion to accretion. The second transition, 
occuring ss 10 6 yrs after the first transition, happens when the torus decreases in density enough to become non-degenerate. At 
this point accretion is no longer effective, and the star's emission is again dominated by vortex expulsion. The left-over non- 
dgenerate disk is slowly consumed resulting in RRAT-like behaviour. Shown in the right panel are the three RRATs studied here, 
two of which we argue have spent long enough time in the tours (accretion) era before transiting back to having their emission 
dominated by vortex expulsion. The arrows depicts the large change in their X-ray luminosity during the transition. Also, in the 
right panel are the 2 knowns AXP transients (AXP IE 1547-54 and AXP XTE J18 10-197), which fall within the first transition, 
from vortex expulsion to accretion dominated emission. 



Rm = R\ + A/? t at Ttorus- We can estimate AR t by integrating 
dt/dr = r B /Ar w up to r torus , with r B and Ar w given by eq(5) 
and (17) in OLNII, respectively. We get 



95 



R 



205/24 
L15 



* 1 + 



ARtj 



73/24 



which leads to 2 cases 
31 



15 - 



R 



ARt] 



181/24 
t,15 

31 



- 1 



IfR U 5 > 1.7 



R. 



132/73 
t,15 



IfRti5 < 1.7 



(A.l) 



and viscous spreading), until it reaches non-degenercy densi- 
ties, at which point accretion becomes much less efficient, and 
the sources luminosity becomes dominated by vortex (mag- 
netic flux) expulsion from the star. This results in a decrease 
in luminosity, as is illustrated in figure |2]by the vertical lines. 
It is this transition, we argue, that is responsible for producing 
RRAT behaviour. 



(A.2) B.2. Transition from vortex to accretion dominated era 

Before the torus phase, there exists a critical P at which a transi- 
tion from vortex dominated to accretion dominated luminosity 
will occur. This critical value is found when L v L acc , 



Appendix B: The L x -P diagram in our model 

B.1. Transition from accretion dominated to RRAT era 

The left panel in Figure |2]illustrates the torus and RRAT eras in 
the Lx-P diagram. In the torus phase (lasting roughly 10 5 yrs; 
eq. [5]) the X-ray luminosity is dominated by accretion from the 
torus, during which time the source spins down due to mag- 
netic braking, thus, evolving horizontally towards smaller P. 
As the torus is accreted its density decreases (due to accretion 



SxKr^ss" 1 



''X,0.1 QS.1.4 



(B.l) 



which implies 10 12 < .P cr i t (s s l ) < 3x10 11 for typical values 
of the torus radius in our model 1 1 < R t (km) < 32. 

The two known transient AXPs (AXP IE 1547-54 and 
AXP XTE Jl 8 10- 197) fall within this transition region (see 
right panel of fig.|2]i. In an upcoming paper we investigate tran- 
sient radio-emission mechanism associated with these sources. 
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B.3. Statistical appearance of SGRs 

In our model, features and statistics of SGRs can be accounted 
for if they are born (via the quark-nova mechanism) from 200 
ms neutron star progenitors. The statistical lin^| for 200 ms 
(shell harboring) sources in our model is, with a birth rate of 
one per 500 yrs (as estimated by Leahy & Ouyed (2007b) I, 



fsmt.,200 ~ 1-2 S fi^l 3 5 P a200 ms^QS,10 M Qs!L4 ( R2 ) 
ntat.,200 ~ O X 1U SS « ,15^0,200 ms K QS,10 M QS,1.4 • 

If we adopt an initial magnetic field strength of Bo = 10 16 G 
for the 200 ms objects we find a statistical line at P s m.,200 ~ 5.3 
s and Pstat.,200 ~ 2.7 x 10~ 10 s s _1 , which is to the right of 
SGRs location in the Lx-P diagram. In fact for these objects, as 
they evolve to lower P the shell is consumed and they becomes 
less active. Thus our earlier suggestion (OLNI and OLNII) 
that what is referred to as AXPs in the literature are old shell- 
harbouring sources in our model that consumed most of their 
shell. 



B.4. Statistical appearance of AXPs 

The statistical line for millisecond (torus harboring) sources is 
given by eqs. (0 and ([8]) in our model (for t ~ 500 yrs): 

P staL)1 ~ 0.2 s CCWu ( B - 3 ) 
P , ~ lO-'W R 2/i P 1/3 R 4/i M~ 1/3 

However as can be seen in the right panel of Figure|2] all AXPs 
are to the right of the above derived statistical line. As we 
will argue below, all of the observed sources to the right of 
the 10" s s _1 line are SGRs (shell harboring quark stars in 
our model) with different shell masses with older ones having 
smaller shells thus with fewer and less energetic bursts (see 
OLNI). 

Torus-harbouring sources in our model are first seen at pe- 
riods of ~ 0.2 s and period derivative of the order of 10~' 1 s s _1 
We thus predict a population of AXPs with sub-second period 
and an X-ray luminosity ~ 10 34 erg s -1 , from eq(Q~|i. 



B.5. SN transparency line 

For completeness, we also define a transparency line related to 
the time needed for a supernova remnant to be transparent to 
X-ray. A rough estimate gives 

?SN,transp. ~ 2 yr MsN,l()Vsjy[,4000 ' (B-4) 

for a SN ejecta of mass 10M G and velocity of 4000 km s~'. 
In our model, for AXPs (fo.ms ~ 1 and Z?o,i5 ~ 1) we get 
^SN.transp. ~ 39 ms and PsN.transp. ~ 2.3 x 10~ 10 s s" 1 . For SGRs 
(-Po.ms ~ 200 and B ,i5 ~ 10), we find -P S N,tran Sp . ~ 1.1 s and 

^SN.tiansp. ~ 6 X 10~ 9 S S" 1 . 



2 The statistical line is the age of the youngest observed source in a 
population: age ~ 1 IB for birthrate B. 



